function res = fn_G(x, T)

if x/T>33
    
    res = 0;
    
else

    res = 1./(1 + exp(x/T));
    
end

end
